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Abstract 


This document is a final technical report on a National 
Aeronautics and Space Administration research gtant titled 
”The Physics of Turbulent Plow", NSG 3280, June 1, 1979 to 
May 31, 1980, from NASA Lewis Research Center. The investi- 
gators proposed to use shock tube discontinuities so as to 
test and extend a kinetic theory of turbulence. In shock 
wave and contact surface fluctuations, coherent phenomena 
have been found which provide new support for the microscopic 
nonempirical approach to turbulent systems, especially those 
with boundary layer-like instabilities. 
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I . Introduction 

The need to develop a complete theory of turbulence has 
become urgent. The current potpourri of empirical data fits, 
phenomenological models, and computer based simulations are 
becoming more and more unacceptable due to escalating costs 
and increasingly inadequate achievements. More importantly, 
they share a fundamental deficiency: traditional approaches 

violate form invariance and are ultimately inconsistent with 
the basic laws of physics. 

Recently, by contrast, a kinetic theory of turbulence 
derived from basic principles and without empirical constants 
has been developed. In this theory, with the assumption of rigid 
molecular chaos removed, turbulent effects become an explicit 
part of the fluid equations; expressions are obtained for 
familiar quantities such as fluctuations in spontaneous stress 
and heat flux, Reynold's stress and turbulent heat flux. 

Other effects such as reaction rate distortion in turbulent 
flows are predicted. We have therefore sought to test some of 
these expectations using the boundary layer-like discontinuities 
in shock tube flow. In this report, we review some of our 
activities during the year of NASA support and speculate on 
their implications. 
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II. Density Measurements form Crossed Beams at High Extinctions 

The importance of intersecting beams of radiation as 

an alternative to hot wire anemometry has already bean shown 

by Fisher and Krsuse^. We wish to extend tlieir analyses so 

that the technique might be applied to the unstable contact 

surfaces which we have recently reported , For this, we use 

a dilute mixture of N 20 ^ ^ 2 NOg in N 2 strongly absorbing 
0 

at 4350 A; other experimental details can be found in 
reference 2. 


Referring to Fig, 1, the beam intersection point has 
coordinates (x,y,z) ; n, and x define distances from this 
point in the x,y and z directions. If K is the extinction 
coefficient for a given species in the turbulent flow, then 
K is a function of position and time and the statistical 
fluctuations in K will yield statistical fluctuations in a 


flow property, 
at time t is: 


For beam the intensity recorded at D 


1 


t 

Il(t) = Iq exp{ - K(x,y + n,z, t) dn) 


The instantaneous extinction coefficient written as the sum 


of its time-average mean value and a fluctuation relative 
to this value is: 


K(x,y 2, t) = t k(x,y + n, z, t) . (2) 

Further, determining = Iq exp (-Kj^ ^ ) , eq. (1) becomes 

I^(t) = <I^> exp{-j k(x,y + n , z,t) dri). With similar 
treatments for S 2 D 2 , we obtain: 



( 3 ) 




Figure 1. The cross-beam correlation experiment. The circle repre- 
sents the cross section of our shock tube. The unsteady 
flow is into the plane of the paper. S 2 ^,S 2 /Di.D 2 are the 
two sources and detectors of light (X = 4350 A) from 
the Hg-vapour lamp. The output of the photomultiplier 
tubes (PMT) is fed into the A.D.C. (BIOMATION 8100) 
and then to the PDP-10 via interface. 




►31 h- 


4 . 


Eqs. (3) are the modified detector currenta which are fluctu- 
ations about a local average level specified by <l2> and <l2^* 
Taking the logarithmic product of egs. (3) and integrating 


over T 


fT- 

in • in i2<t) dt 



Jc(n,y+n# 


z rt) • 


k(x,y, z +T,t)dt} 


- F(x,y,z) 


( 4 ) 


F(xty,z) is the covariance of the fluctuations along and 

S2D2* Assuming a stochastic nature for the fluctuations > the 
important contributions to F(x,y^z) will come from a correlated 
area about the beam intersection point. Defining this 
correlated area as we obtain: 

l*T 2 

P(x,y,z) ■ ini^(t)»An X2(t) - k (x,y,z) A . ( 5 } 


The new relationship above is similar to the one derived 

in reference 1; however^ we have avoided their 

restricted forms for i2(t) and i2(t) which would only apply 

to weak extinctions. The measured quantity F(Xfy,z) is 

directly proportional to the local turbulent intensity at 

the point of beam intersection. 

Using this approach; the gross features which we observe 

in the point density fluctuations in the NO2 driver gas 

mixtures are the same as those previously reported from 

pressure; interferometric; and line-averaged absorption 
2 

measurements . Bursting behavior is clearly indicated in 
the k (X;y;z) histories, in addition, power spectra for 
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the fluctuations all show a dominant mode at about 400 Hz 
with increasing harmonics at increasing flow velocity. 

These results are illustrated by the data in Figures % and 3. 
The emergence of a unique "characteristic" frequency tor 
unstable contact surfaces is reminiscent of similar effects ' 
in turbulent bursts associated with wall boundary layers 
and free jets at transition^ ' ^ . This reinforces the sense 
of phenomenological kinship between the contact surface 
and other discontinuities for which the boundary layer 
approximation seems appropriate^. However, the underlying 
physical significance of these frequencies is not yet 
determined. 
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Figure 2. Crossed Beam Data from the Turbulent Contact Surface, 

(a) Samples of Input Data. The first and thir^ traces are 
the absorption of beams S^jD^^ and S2D2 at 4350 A. The arrival 
of the contact surface with NO2 occurs about 450 ys after 
the display trigger signal. The second and fourth traces are 
the absorption at S^Di and S2D2 normalized to the secular 
trends in the first and third traces respectively. Firing 
conditions were^ : the driver gas mixture was 2% N2O4+2NO2 

in N2 at room conditions; the driven gas was N2; the shock 
wave Mach number, Mg, was Mg = 1.6; and the gas velocity at 
the contact surface was U2 * 313 m/s. The signals shown 
are 1240 ys long in laboratory time; the vertical scales 
are arbitrary (linear) units. (b) Samples of the Point 
Density Fluctuations. (See text, equation (5).) The upper 
trace is the crossed beaun covariance for the data in (a) . 

The lower trace is that for a similar firing at Mg * 3.2 
for which the contact surface arrived at a time roughly 
coincident with the display trigger. 
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Figure 3. Turbulent Power Spectra at the Contact Surface. The upper 
and lower traces are from the corresponding data in 
Figure 2 (b) . 
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III. Turbulence in Argon shock Waves 

Gasdynamic discontinuities are known to be inherently 
unstable when the pressure gradient and the density gradient 
are opposed^. For example, detonation waves show turbulent 
bursts as a result of the post-ignition expansion of the 
burning mixture of gases , Contact surfaces provide turbulent 
intermixing of the otherwise separated gases when a large 
increase in density associates with a perturbative decrease 
in pressure at the front . In these examples, the vorticity 
required for turbulence comes from curvature at the discontin- 
uity^, An ionizing collisional discontinuity can also show 
5 

curvature ; a variety of mechanisms can be suggested which 

6 7 

might cause the growth” and enhanced manifestation of random 
fluctuations in it. Therefore, we have chosen to look for 
turbulent effects in ionizing shock waves in argon. 

Our shock waves are produced in a cylindrical arc dis- 
charge stainless steel shock tube of 5 cm diameter and 185 cm 
length. The tube is pumped down to less than 10 Torr before 
filling. The shock tube is filled with high purity room 
temperature argon and fired without the use of a diaphragm. 

No use is made of any external magnetic fields. The discharge 
is generated by triggering a 14.5 yP capacitor charged to 
18,5 kV and typically takes about 50 psec for completion. The 
shock wave's speed can be varied between 5 x 10 cm/sec and 
3 X 10^ cm/sec by varying the initial filling pressure between 
250 mTorr and 25 mTorr. The test section is 160 cm from the 
discharge. Two piezoelectric pressure transducers with 15 cm 
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separation are positioned fore and aft of the test section's 
window. At the window, a pair of ienses and a pinhole diaphragm 
are used to collect white light emission (with a 3,5 mm x 1 mm 
spatial resolution) , At the same location, two steel probes 
of 0,75 mm diameter are placed in the tube with varying separa- 
tion along the axis of the tube. The probes are biased at 
negative voltage so as to collect ion current. 

The pressure transducer signals and the white light 
emission signals are used to determine the velocities of the 
shock waves which we produce. The white light signals are 
also used to confirm the general qualitative features of the 
data from the various probe signals. The probe signals are 
sent simultaneously to a transient ADC recorder (Biomation 
8100) which has a 25 MHz bandwidth in the chopped mode. The 
digitized data from the probe current are used to obtain the 
plasma power spectrum via a Fast Fourier Transform, The phases 
of the Fourier transforms are compared; the wave number, k, 
is obtained using Aa =: ik*(Ay) = ik_Az where Aa is the phase 
difference and Az is the separation of the two probes . 

With these several diagnostics, we have confirmed the 
existence of turbulent-like fluctuations in collisional shock 
waves . Figure 1 illustrates this with the probe signals from 
three different shock tube firings under three different 
conditions. The arrival of the shock wave is observed as a 
sharp rise in the probe current. In each case, rapid irregular 
oscillations in the signals are found which persist for 5 to 
25 v>secj they are then followed by a relatively quiet decay 
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Figure 1. Sample Data: 

Probe Signals. 

Top 

Trace : 

#051606A, 

P = 

100 mTorr, M = 

16. 

Middle Trace: 

#051503A, 

P = 

60 mTorr, M = 

26. 

Lower Trace: 

#043010A, 

P = 

40 mTorr, M = 

36. 

(M is the Mach number 


of the shock wave and p is the pressure in the shock 
tube prior to the discharge of the capacitor.) Labor- 
atory time, tjj, and distance, dg, plotted above are 
related by dg^t^ x V_ where Vg is the flow velocity 
at the shock front. ^ ^ 
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phase. The data are shown in a frame of reference fixed with 
respect to the gas at the shock wave. The general features 
of the fluctuations appear to be Mach number independent; that 
is, the fluctuations are essentially frozen in the moving frame, 
lasting from 12 to 15 cm, with the same overall shape. 

Fourier analyses of the fluctuations in probe signals 
provide further evidence of a frozen pattern. In Figure 2a 
are shown the power spectra of two probe signals (designated 
A and B) with probe B 0.75 cm downstream from probe A along 
the shock tube's axis; the pressure in the shock tube is 30 
mTorr and the shock wave is moving at Mach number M = 42 at 
the test section. The third spectrum, designated B', is 
obtained from a shock tube firing at p = 40 mTorr and M = 36. 
The frequencey axis for B' is adjusted by a factor of 1.16 
in accordance with the change in Mach number. We notice that 
there is a general peak to peak correspondence in A ai/d B. 

This suggests a longitudinal reproducibility in the statistical 
behaviors of the shock front fluctuations. The spectrum B' 
suggests that the fluctuations in two different firings under 
similar conditions are similar in their statistical features. 

In addition, we find that the phase velocity of each mode is 
about the same as the group velocity. This is illustrated in 
Figure 2b for the two firings in Figure 2a; the dispersion 
relations for A and B and for A' (not shown) and B' are nearly 
linear. In both cases, the group velocity is slightly lower 
than the measured velocity of the shock wave (Vg/Vg'Vi 0.9). 



Figure 2. Sample Data; Spectral Analyses. (a) Power Spectra. Upper 
and middle traces are from two probes (A and B) separated 
by 0.75 cm along the axis of the shock tube for a firing 
with p = 30 mTorr and M = 42. The lower trace (B') is from 
a firing with p = 40 mTorr and M = 36. The vertical scales 
are in arbitrary (logarithimic) units. The horizontal 
scale is frequency (MHz) . (b) Dispersion Relations. The 
lower trace is for signals A-B above. The top trace is for 
signals A'-B' for which B' is shown above. 
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However, v is approximately the same as the local flow speed 

y 

behind the shock front. Therefore, in the frame moving with 
the flow velocity, the group velocity is essentially zero,- 
as required for a frozen fluctuation pattern. 

The spatial evolution (in the moving gas frame) of 
selected modes in the probe signals has been obtained by the 
inverse descrete Fourier Transform after filtering by a 
"brickwall” type, zero phase, digital bandpass filter with a 
variable bandwidth. Samples of these data are shown in Figure 
3. We find that the spectra are quantized and that the funda- 
mental mode has a wave length of about 2.5 cm, i.e., about 
one-half the diameter of the shock tube. All the modes are 
heavily damped; however, the higher the frequency, the faster 
the growth and damping rates. Finally, we notice that the 
envelopes of different modes oscillate with different periods, 
thereby suggesting the. possiblility of intermodal energy 
transfer. 

From these measurements , we see that the fluctuations in 
our shock waves display typical turbulent signatures. The 
presence of a dominant mode in the fluctuation spectra which 

0 

we find is appropriate for turbulence in its earliest phases . 
The reproducibility of the characterizing modes in the fluctu- 
ations which we have shown is to be expected for gas systems 

. 9 

where the boundary conditions- are unchanged . The use of a 
point diagnostic provides us evidence for turbulence in a 
shock front which is more conclusive than that available from 
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Figure 3. The Spatial Evolution of Selected Modes: Run 

#051606 A. Top trace: fundamental mode, f jvO.18 MHz; 

0.15if<0.22 MHz; Middle Trace: second harmonic, 

f2'^’Q.34 MHZ; 0,26if^0.41 MHz. Bottom Trace: fourth 

harmonic, f^'^^0.72 MHz; 0.62s:f^0.83 MHz. The location 
of the leading edge of the shock wave as seen in the 
raw probe is indicated by an arrow. Horizontal axis: 
distance (cm). Vertical axis: Amplitude squared 

(arbitrary units) . 
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line-averaged techniques^®. The analyses of fluctuations v/hich 

we report avoid the bandwidth limitations of earlier reports on 

nonuniformities in shock waves which used hot-wire anemometry^^. 

We also have shov;n that the recently proposed argon shock front 

12 

morphology which ignores fluctuations is incomplete 
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IV. Unsteady Turbulent Shear Flow in Shock Tube Discontinuities 
In collisional fluids and plasmas, there are several 
mechanisms through which perturbations in unstable flow may be 
amplified: (1) Rayleigh-Taylor instabilities, arising when a 

pressure gradient is opposed by a density gradient? (2) rotationa*: 
convection instabilities, arising when vortices are produced in 
such a manner as to enhance the spatial dimensionality of the 
flow? (3) ion acoustic instabilities , arising from the direct 
coupling of the induced electromagnetic field of the perturbation 
and local ion dynamics? (4) reaction-diffusion instabilities 
arising becasue of an inhibition of dispersion by a temporally 
competitive non-equilibrium process. The first two are tra- 
ditional "fluid” mechansims , the second two being usually 
associated with with highly ionized gases. These natural insta- 
bilities can lead to turbulence in boundary layer-like flows 
which show two kinds of statistical fluctuations [1] : (1) random 

fluctuations with indeterminate phase relationships between 
various flow segments? (2) localized fluctuations with deter- 
ministic dynamics and fixed phase relationships from flow segment 
to segment. The second category (i.e., flow with "bursts" or 
"large coherent structures") has recently grown in importance 
since it plays a large role in the transport and mixing processes 
from which turbulence derives much of its crucial significance [2] . 

There is no complete theory of turbulence for these phenomena 
[3,4,5], However, approaches to the problem exist which sug- 
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gest that (near transition) reaction rates can be distorted [6J 
and that a boundary layer approximation affords a generalizable 
quantitative context for measurements on turbulent bursts in 
unstable contact surfaces {71 and detonation waves 18]. Further- 
more, we have just completed a study of bursting instabilities 
in collisional ionizing shock fronts which has analogous impli- 
cations [9J . All of the results just cited have treated the 
instantaneous in-place manifestations of turbulence. As an 
extension of these interests, we report here on the real time 
evolution of two kinds of unstable shock tube discontinuities, 
both of which are boundary layerlike due either to the nature 
of the velocity discontinuity [10) or to the presence of 
curvature [11] . 

The flow properties of the shock tubes which we have used 
are summarized by the x-t plot in Figure 1. Energy is stored in 
a driver section: for the arc discharge tube a 14.5 )Jf capac- 

itator is charged to 18.5 kV; for the pressure ruptured tube, 
a diaphragm contains a pressure load of between 1 to 3 atm. The 
energy is released by the sudden discharge (^20 usee) of the 
capacitor in the arc-discharge shock tube and by the sudden 
rupturing of the diaphragm in the pressure loaded shock tube. 

This release causes a blast wave, or alternatively, a contact 
surface to propagate downstream into the driven section pre- 
ceded in both cases, due to the strength of the disturbance, 
by a shock wave. The details concerning the production of 
shock waves by shock tubes are discussed in many places [12] . 
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In the pressure ruptured tube, dur measurements are performed 
on the contact surface. A firing to firing variation in flow 
velocity at the turbulent contact surface is available through 
the variation in primary shock wave Mach numbers, M_, 

*9 

1.5<Mg<3.0 producing a range in local Reynolds number [71 of 
100<Re<1000. The test section contains two ports separated 
by 15 cm at which crossed- beim stations are placed. Blue 

O 

light (4350 A) absorption of NO 2 is used and a direct measure 
of the covariance of density fluctuations at the point of 
intersection is obtained from [13] : 


F(x,y,z) 



Jin ij^ (t) • in i, (t) dt = I(x,y,x)A 


where ij (t) is the light intensity at path (j) with respect 
to the instantaneous average level <G •> (i • (t) /<G . (t) >) , 

»J J 

T is a convenient averaging interval, A is the cross section 
at the intersection of the beams, and I(x.y.z) is a measure 
of the intensity of turbulent density fluctuations. In the 
arc driven tube, our measurements are performed on the shock 
front. Over the length of the test section, the shock wave 
* is heavily attenuated, providing a range of Mach numbers 

during a firing which can vary by as much as 5%/cm, Conven- 
tional electric probes [14] operated under ion saturation 
conditions are the major diagnostics. The probes are separated 
by 13 . 2 cm and give a direct measure of ion density through- 
out the plasma of interest. In both cases, the primary diagnostics 
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Figure 1. x-t Diagram for Shock Tube Discontinuities. The 
arc-discharge tube (upper sketch) is 190 cm long 
and 5 cm in diameter. The test section is 160 cm 
from the 10 cm long discharge section with inlets 
for electric probes as indicated above. The 12.7 
cm diameter pressure ruptured shock tube (lower 
sketch) has a 153 cm long driver section and a 336 
cm long driven section. The test section is 
153 cm long, immediately adjacent to the diaphragm. 
See references [7] , [10] , and [13] for details. 
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are supplemented by the usual assortment of pressure gauges 
and line-averaged optical sensors. 

First, we turn to the question of the dynamics of bursts 
in these two systems. Direct observations using pressure 
gauges roughly 3 cm apart have already provided evidence 
that turbulent bursts in contact surfaces propagate as dis- 
crete flow entities with velocity at roughly 90% of the local 
flow velocity [7]. Similarily, turbulent bursts in collisional 
ionizing shock fronts measured over distances of about 0.7 
cm also move at velocities of roughly 90% of the local vel- 
ocity; in this case, the velocity measurements are derived 
from the dispersion relations calculated from «he power spectra 
[9] . -These behaviors are generally consistent with the reports 
on large coherent structures in other manifestations of 
boundary layer flow [15] . 

A sample of the contact surface data is given in Figure 2. 
The driver gas is a N20^ mixture {1.0%(N2O2 :+ 2 NO 2 ) ; 99% He). 

In Figure 2 (a) , the arrival of the contact surface is indicated 
by the sudden increase in absorption of blue light at the 
first and second stations. The lag in the two signals results 
from the time of flight between the two stations. Figure 2(a) 
also shows the prominent distortions in what should be other- 
wise a smooth monotonic change in the light signals. These 
are the turbulent bursts as they show up in line-averaged 
blue light absorption. In Figure 2(b), the covariance of each 
station is displayed. The signatures of the turbulent bursts 
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!a) 


(b) 






Figure 2, Crossed Beain Data; (N2O4 ^ 2N02;He) Mixture for 

Driver Gas; N2 for Driven Gas. (a) Light Absorption. 
The upper trace is the absorption in one of the two 
crossed beams at the first station (Station A) . The 
lower trace is the absorption in one of the two 
crossed beams at the second station (Station B) 15 cm 
downstream. Mg a * 2.1; Mg n * 2‘.3. The vertical 
axes are in arSltrary (amplitude) units; the horizontal 
axis is sample number, 4 ys/sample. (b) Normalized 
Covariance I (x,y,z) for Crossed Beams. This the 
measure of turbulent density fluctuations, at the 
points of intersection, for the two stations. The 
vertical axes are arbitrary (power) units; the 
horizontal axes is local particle time in msec. 

(c) Power Spectra for Crossed Beams. The vertical 
axes are arbitrary (power) units. The horizontal 
axis is sample number, 122 Hz (sample) . 


here are rhe resuXts in the region o£ intersection, a volume 
of approximately (0»5 cm) . Power spectra for the covariance 
function are performed using standard FFT digital analyses 
and are displayed for these data in Figure 2(c). These three 
sets all imply a particle-like evolution of bursts. The 
velocity of the contact surface increases by more than 10% 
between stations A and B; however, the power spectra show 
a substantial repeating of the overall features. The results 
represented in Figure 2 are found in all cases studied. 

A set of ion density profiles from the arc-driven shock 
tube is shown in Figure 3. In the top two curves, the shock 
front data from a single firing are given in which the argon 
pressure prior to discharge was 40 mTorr; the shock wave 
decelerates between the two stations (with 13.2 cm separation). 
In both cases, we can identify bursts in the plasma; we can 
also observe that the disturbances persist over the 13.2 cm 
distance. Peak to peak correspondence of the clumps can be 
easily traced as indicated by the assigned labels. This allows 
us to notice that the bursts seem to catch up with the shock 
fronts, a behavior of a sort consistent with a soliton-like 
dynamics [16] . The results represented in Figure 3 are found 
in all our data suggesting a particle like behavior under widely 
varying environments for the evolution of bursts in turbulent 
shock fronts . 

Secondly, we have examined the evolution of the spectral 
behaviors of the turbulent bursts. Figure 2(c) suggests a 
persistance of pratdnent components of the spectral profile 
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for an accelerating discontinuity. In Figure 4, the power 
spectra are shown for a firing in which the contact surface 
decelerates by roughly 20% between the first station and the 
second. The gross features of the spectral profiles in 
terms of prominent modes are certainly reproduced in the 
evolution of the bursts. In firings processed to date, we 
find this to be the case, including the apparent tendency for 
a burst during its history to provide more high frequency modes, 
the higher its velocity. 

For bursts in the shock fronts, the same general features 
are observed. By way of illustration, a firing is displayed 

in Figure 5 in which the shock front decelerates and the 

. . “I 39 

spectral profiles are determined. The behaviors ‘ ' and 

"■1 67 

SgCK ' indicate a turbulent-like spectra dependence [17] ; 

fi 

SaK^. However the differences in power law trends indicate 
that the relative strength of the higher wave number modes de- 
creases as the shock front propagates and decelerates. 

Filially, we turn to the possibility of an explicit con- 
nection between these bursts and other boundary layer flow. 

In this, we can be guided by the well established consensus 
v;ith regard to the appropriateness of the treatment of contact 
surfaces as boundary layers [18j and in shock-tube discontinuity 
[8] with the Orr- Sommerfeld solutions for Falkner-Skan velocity 
profiles . 

In the latter regard, the amplifications of fluctuations 
is expected to follow a predictable serpentine profile, para- 
meterized by a characteristic frequency and a maximum in 
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the turbulent intensity vs Reynolds number plane [191. This 
behavior is found for turbulent bursts in a boundary layer 
[20,181 and for turbulent bursts in a detonation wave [11,81. 
Although our earlier results on contact surfaces are suggestive 
[7] ; the pressure gauges used are wall mounted and band-width 
limited; the interferometer provides line-averaged (rather 
than point) data; and the Mach number range studied does not 
provide adequate determination of the change in turbulent 
intensity with increasing Reynolds number. However, our 
present data are based on point diagnostics, reach higher 
flow speeds with 4 times the previous bandwidth and permit us 
to fashion a test of the predictions of Orr-Sommerfeld-like 
behaviors. Specif icially; (a) the strength of the fluctuations 
in a burst should begin to decrease with increasing Mach number 
Mgf^2 (see figure 12 in reference 7 and associated discussion 
for the relationship between shock wave Mach number and contact 
surface Reynolds number) ; (b) a single dominant frequecny should 
characterize the spectral profiles of fluctuations associated 
with the turbulent bursts. 

Our data for turbulent bursts at the contact surface confirm 
these expectations. The local reynolds number is given by 
Re a £ where £ is a characteristic length scale; £ a tj^ where 
tj^ is the separation between the burst and discontinuity [8]. 

Thus , one can map I vs into I vs Re where I is the turbulent 
intensity. The expectation (a) above means d I/dt < 0. 

Figure 6 (.a) shows a plot of AI /At,, where, in each covariance 

X ij 

with at least two bursts, a computation of (I .-I .)/ 
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Figure 3. Electric Proble Data: Shock Heated Argon. The upper 

two traces give relative ion density during a dis- 
charge into Argon at initial pressure p]_=40 mTorr; 
Mg^;^=45, Mg 3=33/ The lower traces are for P]_ = 30 
mTorr; 52, Mg g=39. The vertical axes (not 

shown) would have be4n arbitrary (amplitude) units; 
the horizontal axis in in cm; measured in a frame of 
reference at rest in the local gas. The labels (a,b, 
c,d) suggest corresponding clumps in the ion density 
profiles. 
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Figure 4. Crossed Beam Power Spectra; (N 204 ;^: 2 N 02 ;N 2 ) Mixture 
for Driven Gas. Mg ^=2.35; Mg g=2.2. The scales 
here are the same as in Figure 2 (c) . 
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Ion Density Profiles and Associated Power Spectra. 
P2^=60 mTorr Argon? Mg ^=33, Mg b= 25. The scales 
for the upper two traces are tAe same as Figure 3 . 
The lower two traces are the corresponding power 
spectra. The vertical axes are arbitrary log- 
arithmic units; the horizontal axis is wavenumber 
K in cm”l from K = w/Vg^s • 


Figure 5. 
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Figure 



(a) (b) 


Crossed Beam Data on Boundary Layer-Like Profiles . 

(a) Ijj. is the value of I at the peak of an identi- 
fiable burst and t^ is the time interval between 
that peak and the arrival of the contact surface 
at the same point in a frame of reference at rest 
with respect to the local flow (see Figure 2) . For 
two peaks in a given covariance profile AIjj/AtL = 
(Ijj, 2 ^x,l) / 2“^L 1^* Each profile corresponds 
to a definite value for Me* All values of AIx/At^ 
obtained at comparable values for Me are averaged 
without regard for flow history. (b) Each identi- 
fiable peak in a power spectrum, f-., is interpreted 
in terms of a relationship where 

nj^=l,2,3... and f^ is the estimate of f^ , the funda 
mental frequency , provided by fm,i» The values of f 
obtained from comparable values of M^ are averaged 
without regard for flow history, 




fi>. 
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t, j-t- .) is performed; averages are determined at a given 
Mg and the profile AI^^/At^ vs Mg is achieved as a representation 
of Al/ARe vs Re. The negative slope is qualitatively consistent 
with the Orr-Sommerfeld profiles. In addition, in each of the 
power spectra, an estimate of the fundamental mode, f^, is 
made by assuming an integral relationship between it and each 
prominent spectral component. This given an estimate of f^ 
at each firing and in all spectra when we require that a, single 
fundamental frequency should describe all data. The results 
are given in Figure 6 (b) . The average value of the fundamental 
frequency is f^ % 1041 Hz; the data are not inconsistent with 
the assumption of a constant value for this frequency over the 
range of speeds studies. 

With these results, we see the first evidence of the 
possibility of nonempirical unified treatment of all boundary 
layer-like unsteady shear flow without regard for flow history. 
In particular, we see new manifestations of a "constant" 
characteristic frequency with these turbulent bursts of the 
sort previously found in conventional boundary-layer flow [21J, . 
The usefulness of the interpretation of turbulent bursts as 
characterized by Orr-Sommerfeld profiles should therefore be 
tested for the widest possible range of flow environments . 
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V. Summary 

Our investigations have produced several new results. 

The formulation of the crossed beam correlation technique 
has been generalized to include strongly absorbing media. 

The first measurements of point density fluctuations at con- 
tact surfaces have been obtained. The presence of turbulent 
bursts is confirmed and a characteristic spectral frequency 
of approximately 400 Hz is estimated. Irregular density 
fluctuations with turbulent-like behaviors are found in ion- 
izing shock fronts produced by our arc driven shock tube. We 
use electric probes as the primary disagnostic. Spectral analyses 
show statistical patterns which seem frozen-in and charactizable 
by a dominant mode and its harmonics. Using a pressure-ruptured 
shock tube and an arc driven shock tube, we have studied the 
evolution of turbulent fluctuations at contact surfaces with 
*^204 t ionizing shock fronts in argon, 

we have focused on point density diarhostics derived from crossed 
light beam correlations and electric probes. Turbulent bursts 
are found for which dynamical and spectral segments with a 
unique and characteristic frequency, independent of flow history 
and overall flow conditions. 
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